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1 Introduction 

At every point, a smooth plane curve can be approximated, to second order, 
by a circle; this circle is called osculating. One may think of the osculating 
circle as passing through three infinitesimally close points of the curve. A 
vertex of the curve is a point at which the osculating circle hyper-osculates: 
it approximates the curve to third order. Equivalently, a vertex is a critical 
point of the curvature function. 

Consider a (necessarily non-closed) curve, free from vertices. The classical 
Tait-Kneser theorem (HUE! (see also [HI ED]), states that the osculating circles 
of the curve are pairwise disjoint, see Figure ^ This theorem is closely related 
to the four vertex theorem of S. Mukhopadhyaya 8J that a plane oval has at 
least 4 vertices (see again P ITU]). 

Figure ^ illustrates the Tait-Kneser theorem: it shows an annulus foliated 
by osculating circles of a curve. 

Remark 1.1 This foliation is not differentiable! Here is a proof. Let / be a 
differentiate function in the annulus, constant on the leaves. We claim that 
/ is constant. Indeed, df vanishes on the tangent vectors to the leaves. The 
curve is tangent to its osculating circle at every point, hence df vanishes on 
the curve as well. Hence / is constant on the curve. But the curve intersects 
all the circles that form the annulus, so / is constant everywhere. 
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Figure 1: A spiral and its nested osculating circles 



Remark 1.2 The Tait-Kneser theorem has an analog in plane Minkowski 
geometry, see [TT] . 

We will prove a number of analogs of the Tait-Kneser theorem; in each 
case, we will obtain a non-differentiable foliation with smooth leaves. 

2 Osculating Taylor polynomials 

Let / be a smooth function of one real variable. Fix n > 1 and let t G R. The 
osculating (Taylor) polynomial g t of degree n of the function / at the point 
t is the polynomial, whose value and the values of whose first n derivatives 
at the point t coincide with those of /: 



The osculating polynomial g t is hyper- osculating if it approximates the func- 
tion / at the point t up to n + 1-st derivative, that is, if / (n+1 * ) (£) = 0. 

Assume that n is even and / <n+1 ' ) (t) ^ on some interval I (possibly, 
infinite). 




(1) 
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Theorem 1 For any distinct a,b E I, the graphs of the osculating polyno- 
mials g a and g^ are disjoint. 1 

Proof. To fix ideas, assume that / (ra+1 - ) (t) > on /. Let a < b and suppose 
that g a (x) = gb{x) for some x G R. It follows from (JTJ) that 

and hence (dg t /dt)(x) > (except for £ = x). It follows that ^t(x) increases, 
as a function of t, therefore g a (x) < gb(x). This is a contradiction. □ 

The same argument proves the following variant of Theorem ^ Let n be 
odd. Assume that f {n+l \t) ^ on an interval /. Consider two points a < b 
from /. 

Theorem 2 The graphs of the osculating polynomials g a and g^ are disjoint 
over the segment [b, oo). 




Figure 2: Osculating quadratic polynomials of the function f(x) = x 3 

Figure 121 shows the graphs of the osculating quadratic polynomials of the 
function f(x) = x 3 and Figure H3 of the osculating cubic polynomials of the 
function f(x) = x A . 

1 We will not be surprised if this result is not new but we did not see it in the literature. 



3 



Figure 3: Osculating cubic polynomials of the function f(x) = x 4 

3 Osculating trigonometric polynomials 

Let / be a 27r-periodic smooth function, that is, a function on the circle 
S* 1 = R/27rZ. Fix n > 1 and let t G S 1 . A trigonometric polynomial of 
degree n 

n 

gt(x) = c + yj(ai cos ix + hi sin ix) 

is the osculating trigonometric polynomial of the function f(x) at the point 
t if its value and the values of its first In derivatives at the point t coincide 
with those of /. 

Remark 3.1 The osculating trigonometric polynomial always exists. Ac- 
tually, the following more general fact is classical (see, e.g., [Ej). Let 
i = 1, . . . , N, be a system of functions on an interval I such that the Wronski 
determinant of this system is nonzero everywhere on I. Then, for any suffi- 
ciently smooth function g on I and any t £ I, there is a linear combination 
of functions f\ that, at to, approximates g up to the derivative of order N — l. 
This boils down to solving the linear system 

g u \t ) = J2^ 3 \t ), 3 = 0,...,N-1 

with unknowns q, which has a solution due to non-zero determinant. The 
solution depends smoothly on £ - 
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In our case, the functions /j are 1, cost, smt, cos2t, sin2t, etc., and 
N = 2n + 1. The Wronskian of these functions is constant, which can be 
seen by differentiating its columns. On the other hand, the functions are 
linearly independent solutions of a N-th order linear differential equation, 
hence the Wronskian is nonzero. 

Geometrically, we consider iV-dimensional projective space and the curve 
[fx : ... : /at : <?]• The osculating hyperplane of this curve at the point 
to approximates the curve with N — 1 derivatives. The equation of this 
hyperplane is g = J2 c ifh an d this gives the desired approximation. 

The osculating trigonometric polynomial g t is hyper- osculating if it ap- 
proximates the function / at the point t up to In + 1-st derivative, that 
is, if f( 2n+1 ^(t) = gf' n+l \t). Trigonometric polynomials of degree n are an- 
nihilated by the differential operator V := d(d 2 + l)(d 2 + 4) . . . (d 2 + n 2 ), 
where d = d/dx. Therefore gt hyper-osculates a function / if and only if 
(!>/)(*) = 0. 

Assume that the osculating trigonometric polynomials of degree n for a 
function / do not hyper-osculate on an interval I C S 1 . 

Theorem 3 For any distinct a,b E I, the graphs of the osculating trigono- 
metric polynomials g a and gt are disjoint. 

Proof. It is not hard to see that the real number gf' n+1 \t) depends con- 
tinuously on t (indeed, the function g t depends continuously on t in the 
C 2n+1 -metric). 

To fix ideas, assume that f^ 2n+1 \t) > gf n+l \t) for all tel. We will 
show that dgt(x)/dt > for alH G / and all x £ S 1 (except t = x), and this 
will imply the statement of the theorem as in the proof of Theorem ^ 

Since g t is an osculating trigonometric polynomial, one has: 

gi J \t) = f U) (t), j = 0,...,2n. (2) 



Differentiate: 



dg, 



00 



a -"(*)+fl? +1) (*) = / (, " +1) (*) > 



and combine with (J2J) to obtain: 

a (J) a (2n) 

^ (t) = 0, j = 0,...,2n-l; ^ (t) + g? n+1 \t) = f^\t). (3) 
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The function dg t /dt is a trigonometric polynomial of degree n. If this 
trigonometric polynomial is not identically zero, then it has no more than 2n 
roots, counting with multiplicities. If dgt/dt = 0, then {dgt/di)^ 2n '{t) = 0, 
and the last equality in (j3J) implies that g t hyper-osculates. Thus dg t /dt is 
not identically zero. 

According to 0, the trigonometric polynomial dg t /dt already has a root 
at the point t of multiplicity In. Hence (dg t /dt)(x) ^ for x ^ t. By the 
assumption made at the beginning of the proof and the last equality in ©, 
we have (dg t /dtY 2n \t) > 0. Hence (dg t /dt)(x) > for x sufficiently close to 
t, and therefore (dg t /dt)(x) > for all x ^ t. □ 



Figure 4: Osculating linear harmonics of the function f(x) = x 3 

Theorem El is illustrated in Figure 0] depicting the graphs of osculating 
linear harmonics c + a cosx + b sin x for the function f(x) = x 3 . 

Remark 3.2 Theorem El extends from trigonometric polynomials to Cheby- 
shev systems of functions; the proof remains the same. 

Remark 3.3 For n = 1, Theorem El implies the Tait-Kneser theorem: it 
suffices to consider the support function of the curve and use the fact that 
the support functions of circles are linear harmonics. 
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4 Osculating conies, cubics and fractional lin- 
ear transformations 



Fix d > 1 and consider the space of algebraic curves of degree d. This space 
has dimension n(d) = d(d + 3)/2. At every point, a smooth plane curve 
7 can be approximated, to order n(d) — 1, by an algebraic curve of degree 
d; this algebraic curve is called the osculating curve. One may think of the 
osculating algebraic curve as passing through n(d) infinitesimally close points 
of 7. A d-extactic point of the curve 7 is a point, at which the osculating 
algebraic curve hyper-osculates: it approximates 7 to order n(<f); see [T]. 

In this section, we extend the Tait-Kneser theorem to osculating conies 
and osculating cubic curves. We assume that the curve 7 is free from extactic 
points. We also assume that the osculating conies and cubic curves along 7 
are non-degenerate. 

Consider a smooth function / with nowhere vanishing derivative. For 
every t £ R, there exists a fractional- linear transformation g t , whose value 
and the value of whose first two derivatives at the point t coincide with those 
of /; this is the osculating fractional-linear transformation. As before, it 
hyper-osculates at the point t if the third derivatives coincide as well. This 
happens if and only if the Schwarzian derivative of / vanishes: 



The graph of a fractional-linear transformation is a hyperbola with vertical 
and horizontal asymptotes (or a straight line); we refer to these graphs as 
the osculating hyperbolas. Assume that the osculating hyperbolas for the 
function / do not hyper-osculate on an interval /. Let 7 be the graph of / 
over I. 

From the projective point of view, all non-degenerate conies are equiva- 
lent; since our results are projectively-invariant, we assume that the osculat- 
ing conies of 7 are ellipses. In the case of cubic curves, we assume that the 
osculating cubics of 7 have two components, an oval and a branch going to 
infinity, and that the ovals, not the infinite branches, osculate 7. With these 
assumptions, we have the next theorem. 

Theorem 4 1) The osculating ellipses along 7 are pairwise disjoint; 

2) the ovals of the osculating cubic curves along 7 are pairwise disjoint; 

3) the osculating hyperbolas of '7 are pairwise disjoint. 
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Theorem HI case 2) is illustrated by Figure El 




Figure 5: Osculating cubic curves of a spiral 



Proof of Theorem 01 We will argue about cubic curves, indicating the 
difference with the case of conies and hyperbolas, when necessary. 

Give the curve 7 a smooth parameterization, 7(3). Let T s be the oval 
of the osculating cubic curve at the point 7(3), and f s (x,y) = its cubic 
equation. It suffices to prove that the curves T a and r& are nested for distinct 
parameter values a and b, sufficiently close to each other. 

Give the ovals T s a smooth parameterization, r s (t), such that the tan- 
gency point with the curve 7 corresponds to t = 0, that is, r s (0) = 7(s). 



8 



Let F be the map (s,t) i— > T s (t). We claim that, for t ^ 0, this map is an 
immersion. This claim implies that T a and Tf, are nested for sufficiently close 
a and 6. 

To prove the claim, we need the following lemma. 

Lemma 4.1 Suppose that F(s,t) = (x,y). The Jacobian of F vanishes at 
point (s, t) if and only if 



Proof of Lemma. The covector df s is nowhere zero since the curve T s 
is non-degenerate. This covector vanishes on dF/dt, the tangent vector to 
the curve T s . Therefore the Jacobian of F vanishes exactly when df s also 
vanishes on the vector dF/ds. 

Differentiate the equation f s (F) = with respect to s: 



Thus df s vanishes on the vector dF/ds if and only if f s = and df s /ds = 0. 
□ 

Now we need to prove that the system of equations (j3J) has no solutions 
for t ^ and point (x,y) on the oval T s . Both equations in (j3J) are cubic, 
and they are not proportional since 7(3) is not an extactic point (in the 
case of osculating conies, the two equations are quadratic). By the Bezout 
theorem, the number of solutions is at most 9 (and 4, for conies). In the case 
of hyperbolas, f s (x,y) = (x — a)(y — b) — c where a, b and c depend on s; 
hence df s /ds = is a linear equation in x and y, and system (JJJ) has at most 
2 solutions. 

For any parameter value s, the point 7(3) is a multiple solution of system 
(H|). Since the curve {f s = 0} is the osculating curve of degree d for the curve 
7 at the point 7(5), the function s' 1— > / s (7(s')) has zero of order n(d) at 
point s' = s. We can view / s (7(s')) as a smooth function of two variables s 
and s'. This function vanishes on the line s' = s. According to a version of 
the preparation theorem for differentiable functions jHE] (see also [7]), there 
exists a smooth function <p of two variables such that 



dfs_ 
ds 



0, f.(x,y)=0. 



(4) 




fs(7(s')) = (s-s') m <P(s,s') 
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and (f>(s, s) ^ locally near a given value of s. Restricting this equation to a 
line s = const, we obtain m = n(d). Differentiating with respect to s, we see 
that 7^(7(5')) starts with terms of order n(d) — 1 in s — s'. Then ^-(T s (t)) 
vanishes for t — with order n(d) — 1, because T s approximates 7 up to 
order n(d) at 7(5). Hence the multiplicity of the solution 7(3) of system (jlj) 
is n(d) — 1. 

For d = 2 (the case of osculating ellipses), this multiplicity is 4, and hence 
there are no other solutions. For d = 3 (the case of osculating cubics), the 
multiplicity is 8, and there may be one other solution. However, the number 
of intersection points of an oval with any curve is even, and therefore the 
9-th point (if it exists) lies on the other branch. Therefore system has no 
solutions for t 7^ 0. 

Finally, in the case of hyperbolas, the multiplicity of the solution of system 
(@J) at the point 7(5) is 2, therefore there are no other solutions again. This 
completes the proof. □ 

Remark 4.2 It is interesting to compare Theorem 0] with three results on 
the existence of "vertices": a plane oval has at least six sextactic (i.e., 2- 
extactic) points |8. ; a closed plane curve, sufficiently close to an oval of 
a cubic curve, has at least ten 3-extactic points pQ; and the Schwarzian 
derivative of a diffeomorphism of RP 1 has at least four zeros [2j (see also 
[3 ED])- Likewise, Theorem |H] is related to the existence of "flexes", i.e., 
hyperosculating trigonometric polynomials [T4*j . 

Remark 4.3 In fact, the osculating hyperbolas are the osculating circles in 
Lorentz metric [21 . 

Remark 4.4 Theorem 0] does not generalize to osculating quartics. This 
can be seen on Figure |Hl where several osculating quartics for the curve 
x 2 /3 _|_ ^2/3 _ j are ,-[ rawn _ Each quartic in the picture splits into two ovals, 
one being below and one above the curve. One can see that nearby ovals 
below the curve intersect. 



5 Infinitesimal intersection indices 

In this section, we give some more general results that may highlight the 
proof of Theorem 
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Figure 6: Osculating quartics of the curve x 2//3 + y 2 ^ 3 = 1. 



Consider a smooth map F of a region in R 2 to a region in R 2 . The map 
F gives rise to a family of curves. Namely, for any s G R, we have the 
parameterized curve T s : 1 i— ► -F(s, t), where the parameter t runs through all 
real numbers such that (s, t) is in the domain of F. Suppose that the curve 
T s is given locally by an equation f s = 0, which depends smoothly on s. We 
will assume that df s never vanishes (e.g., if f s are polynomials, then we are 
talking about nonsingular algebraic curves r s ). 

Let (x, y) be a point T s (t) on a curve T s so that F(s,t) = (x,y). Define 
the infinitesimal intersection multiplicity of T s at point (x, y) as the order of 
vanishing of the function 

1 1— > Jacobian[F](s,t) 

at point t. In particular, if the infinitesimal intersection multiplicity is zero, 
then the family F looks like a foliation locally near the point (x, y) and 
for parameter values near s (however, the curves from the family F corre- 
sponding to far-away parameter values may also pass through (x,y)). The 
infinitesimal intersection index of a curve T s (in the family F) is the sum 
of local intersection multiplicities at all points of this curve. The following 
theorem is an infinitesimal version of the classical Bezout theorem: 

Theorem 5 Suppose that all curves T s are algebraic of degree d. Then the 
infinitesimal intersection index of each curve T s is at most d 2 . 
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The proof of this theorem is based on the following lemma: 



Lemma 5.1 The infinitesimal intersection multiplicity ofT s at a point (x, y) 
is equal to the intersection multiplicity of the curves T s = {f s = 0} and 
{^gj = 0} at the same point. 

This is a direct generalization of Lemma 4.1. 

Proof. The Jacobian of F is, by definition, det(^j, ^). We have df s (^) = 
0, and neither the 1-form df s nor the vector field ^ ever vanish. It follows 
that the Jacobian of F is df s (^) times a nowhere vanishing differentiable 
function. In particular, the order of vanishing of the Jacobian coincides with 
the order of vanishing of the function df s (^) at the same point, and, by the 
equality 



considered to a curve T s and express them it terms of the local parameter 
t. Then the order of vanishing of the function is, by definition, the 
intersection multiplicity of the curves T s and {ff = 0}. □ 
Theorem El now follows. 

The following statement provides a description of families F that consist 
of osculating algebraic curves to a given plane curve: 

Theorem 6 Under the assumptions of Theorem^ suppose also that there 
is a smooth plane curve 7 parameterized by s and such that 7(5) G T s for 
each s, and each curve T s has infinitesimal intersection multiplicity n(d) — 1 
at the point 7(5). Then T s are osculating algebraic curves of degree d for the 
curve 7. 

This theorem generalizes the well-known algorithm of finding the envelope 
of a family of lines: the envelope coincides with the locus of points, where 
two infinitesimally close lines intersect. 

Proof. Since ^ = on 7, the curve 7 is the envelope of curves T s (this 
follows from the classical description of the envelope). 




with the order of vanishing of -jf-°F at the same point. Restrict all functions 
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Then the function s' t— > f s ( , y(s')) has a multiple zero at point s' = s. By 
the preparation theorem for differentiable functions [HllHEj, we have 

f s ( 1 ( s >)) = (s-sT<P(s,s>), 

where m > 1 is an integer and is a smooth function of two variables such 
that <fi(s, s) 7^ locally near a given value of s. In particular, the curves T s 
approximate the curve 7 up to order m for s in the chosen neighborhood. 
Reparameterize curves T s to make T s (t) coincide with 7(s) for t = 0. Then 
^-(T s (t)) vanishes at point t = with order m — 1. On the other hand, the 
order of vanishing is n(d) — 1, hence m = n(d). □ 

EndNote. After this paper was completed, we learned that E. Ghys has 
been extensively lecturing on the same subject in the recent years. We are 
pleased to acknowledge his priority in this matter; some of his unpublished 
results generalize those given in this paper. E. Ghys and the authors plan to 
further develop these ideas and to publish a more complete and more general 
account. 

Acknowledgments. We are grateful to Dan Genin for discussions and 
for making figures for this paper in Mathematica. Many thanks to E. Ghys 
for historical comments, in particular, for attracting out attention to Tait's 
paper [TH] . 
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